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Abstract 

The forward-backward long-range multiplicity correlations in proton-proton 
collisions are investigated in the model with two independent sources of par- 
ticles: one left- and one right-moving wounded nucleon. A good agreement 
with the UA5 collaboration proton-antiproton data at the cm. energy of 200 
GeV is observed. For comparison the model with only one source of particles 
is also discussed. 
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1 Introduction 

Recently, the pseudorapidity particle density from a wounded nucleon^ was deter- 
mined by analysing the PHOBOS data [2] on deuteron-gold collisions at ^/s = 200 
GeV in the framework of the wounded nucleon [3] and the wounded quark-diquark 
[4] models. The obtained fragmentation function^ has two characteristic features. 
It is peaked in the forward direction and it substantially feeds into the opposite 
hemisphere, as shown in Fig. 1. In Ref. [6] very similar shape of the contribution 

*e-niail: Adani.Bzdak@ifj.edu.pl 

"'^The wounded nucleon is the one which underwent at least one inelastic collision [1]. 

^In this picture all soft particles are produced independently from left- and right-moving 
wounded nucleons. It is very similar to the assumption of independent hadronization of strings in 
the dual parton model [5]. 
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from a wounded nucleon was found at SPS energy of \fs = 17.3 GeV. The possible 
explanation of the main features of the wounded nucleon fragmentation function 
was proposed in Ref. [7] in the model based on the bremsstrahlung mechanism [8]. 
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Figure 1: The pseudorapidity density of produced charged particles from a right- 
moving wounded nucleon at y/s — 200 GeV cm. energy. pR{ri) substantially feeds 
into the opposite hemisphere which implies specific long-range forward-backward 
multiplicity correlations. 

It is interesting to notice that the picture in which the wounded nucleon populates 
particles into the opposite hemisphere implies specific long-range forward-backward 
multiplicity correlations. This problem will be investigated here in the context of 
the UA5 pp forward-backward multiplicity correlation data at ^/s = 200 GeV [9]. 
Namely, we will test the model with two independent sources of particles, see Fig. 2, 
with the wounded nucleon fragmentation function shown in Fig. 1. For comparison, 
we will also study the model in which particles are produced from only one source 
of particles e.g., a single string spanned between two wounded nucleons. 

Our main conclusion is that the model with two independent sources of parti- 
cles and the wounded nucleon fragmentation function extracted from the PHOBOS 

d — Au data is fully consistent with the UA5 pp forward-backward multiplicity cor- 
relation data at ^/s = 200 GeV. At the same time we conclude that the model with 
only one source of particles is in a very clear disagreement with the data. 

Let us emphasize here that the idea of the multicomponent model of soft particle 
production is not new and was successfully applied to the forward-backward multi- 
plicity correlations data by many authors [5, 10-15]. Our model with two sources of 
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Figure 2: The model of soft particle production in proton- proton collisions with 
two independent sources of particles: left- and right-moving wounded nucleons. 
The arrows indicate that each wounded nucleon may populate particles into both 
pseudorapidity intervals with appropriate probabilities. 

particles has very much in common with the dual parton model [5] or the two-chain 
dual model with non-zero asymmetry of each chain [10, 11]. For other approaches 
seeRefs. [16,17]. 

In the next section basic formulae are introduced. In section 3 we present our 
model in detail and derive for pp collisions the analytical expressions for the correla- 
tion coefficient and the functional relation between the average number of particles 
in the backward interval at a given number of particles in the forward one. We also 
discuss the limit of one-source model. In section 4 our results are tested using UA5 
pp forward-backward multiplicity correlation data. Our conclusions are listed in the 
last section, where also some comments are included. 



2 General formulae 

It is convenient to construct the generating function 

H (ZB, ZF) = y] P (ns, np) 4^4^ (1) 

where P (n^, np) is the probability in pp collisions to find Ub particles in B interval 
and Up particles in F interval, see Fig. 2. It is worth to notice that the generating 
function (1) contains all information about the multiplicities in B and F. 

The correlation coefficient (or correlation strength) is defined as 

_ {riBrip) - (ub) {np) 

/ 2\ / \2 ' V^^ 

where ub and np are event by event particle multiplicities in B and F intervals, 
respectively. If the number of particles in B interval does not dependent on the 
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number of particles in F i.e., {nBUp) = {ub) {np) we have 6 = 0. On the other 
hand, if = n^r in every event then 6 = 1 (maximum correlation). Using definition 
(1) the correlation coefficient b can be expressed by the appropriate derivatives of 
the generating function 

d^H dH dH' 



{nsnp) - {ub) {np) 



OzbOzf Ozb dzp 



d^H dH 



dzl 



+ 



dzi 



dH 
dzp 



- ZB=l,Zp = l 

2" 



(3) 



ZS=1,2F=1 



It is also interesting to study the functional relation between the average number 
of particles (ub) in B interval under the condition oi np particles in F interval 

Y^ub^bP {nB,np) 



{ub) 



T.nB^i'^B^np) 



(4) 



where the numerator and denominator can be expressed by the derivatives of the 
generating function (1) 

1 d^^H{zB,ZF) 

r {'iiB, np) = - 

'riB 



y] P{nB,np) 



2S = l,Zf =0 



E 



^, . 19 d^-H{zB,Zp) 

ubP {nB,np) = — 

"s npl ozb 



dz^F 



(5) 

2b=1,z^=0 

In the next section we calculate (2) and (4) in two models of particle production. 



3 Model 

The schematic view of our model is presented in Fig. 2. We assume that in pp 
collisions all soft particles are produced from two independent wounded nucleons^, 
which populate particles according to the fragmentation function"^ presented in Fig. 
1. Additionally, we assume that in pp collisions the multiplicity distribution in the 
combined interval B -\- F is described be the negative binomial (NB) distribution 

^"^'"■^■''> = r(nVi)rw (l) On) • 

where n is the average multiplicity iylB+F and 1/k measures deviation from Poisson 
distribution. It is obvious that n can be calculated as 

^ = / [pR iv) + PL iv)] dr] = 2 PR (77) dr], (7) 

J B+F J B+F 



^The detailed discussion of this assumption and its successful applications can be found in Refs. 

[3,4,6,18]. 

""^In the cm. frame a contribution from the left-moving wounded nucleon pl{v) = Pr{~v)- 
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where pR [r]) and Pl{i]) = PR {~v) ^^e pscudorapidity densities of produced 
particles from the right- and left-moving wounded nucleons, respectively. 

Recently, wc have shown [18] that the generating function (1) in the framework 
of the above-mentioned model may be written as 

H{zb,Zf) = <yl + j lPLBi'^- Zb) +Plf{1- ZF)]j X 

( n 1 ~^/'^ 

{l + j:[pRB{l-ZB)+PRF{l-ZF)]j , (8) 



X 



where pRp is the probability that a particle originating from the right-moving 
wounded nucleon goes to F interval rather than to B (and analogous for prb,Plb 
and Plf), see Fig. 2. These probabilities satisfy the following conditions 

PLB+PLF^i, PRB+PRF^i- (9) 

These numbers can be easily calculated using the wounded nucleon fragmentation 
function. For instance, prf has the form 

If PR iv)dV .-in^ 
Prf = r-r-i-- 10) 

Taking (2), (3) and (8) into account and performing elementary calculations, the 
following expression for the correlation coefficient in the model with two independent 
sources of particles is obtained 

^ _ n{pLBPLF + PRBPRf) /^^v 

niplf + p\f) + Kplf + Prf) ' 

Assuming that intervals B and F are separated enough so that F can be populated 
only by the right-moving nucleon and B only by the left-moving one i.e., plb = 
Prf = 1 and p^F = Prb = we obtain 6 = 0. Thus, we immediately predict 
the noticeable suppression of the correlation coefficient h with increasing distance 
between B and F intervals. 

In the model with two independent sources of particles the relation between the 
average number of particles (jib) in the backward interval S at a given number of 
particles n^? in the forward interval F has the form [see (4), (5) and (8)] 

'<^^l\--- 2{k + npLF) 2Fi{k/2,-nF,l-nF-k/2,0 
npRB{nF + k/2) 2Fi{k/2, -uf, -uf - k/2,^) 

k + npRF 2Fi{k/2,-nFA-nF-k/2,^)' ^ ^ 



where 

^ PLF{k + npRF) 



PRF{k + npLF)' 



(13) 
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and the hypergeometric function 2-^1(0, —np, ^ — c,C) is defined as 

For comparison we also derive the appropriate formulae in the model with only 
one source of particles e.g., a single string spanned between two wounded nucleons. 
These expressions can be easily obtained from (11) and (12) by deactivating one of 
the sources e.g., the left one. In this case p^B = Plf = (thus ^ = 0) and Prb = Pn 
and prf = Pf where pb + Pf = 1- Finally 

^ (15) 



k + pfu ' 
and 



(riB) Up = , - ( ^ + ) • (16) 

k + ppn \ 2 ' 



This closes the theoretical discussion of the problem. 



4 Results 

In the present section we test our results using the UA5 pp forward-backward mul- 
tiplicity correlation data at ^/s = 200 GeV. The measurement was performed in the 
pseudorapidity range of jr^l < 4 for various symmetric (around 1] = 0) forward and 
backward intervals. In this case, taking the model with two independent sources, 
we have p^p = plb = P and prb = Plf = ^ — P, where probability p is calculated 
from Eq. (10). In the model with only one source we always have ps = Pf = 1/2- 

In Figs. 3 and 4 the correlation coefficient h for various symmetric pseudora- 
pidity intervals is presented. The experimental data (squares) are taken from Refs. 
[9,19]. The grey and dashed bands represent the results of the model with two 
independent sources and the model with a single source, respectively. The widths 

of the bands reflect the uncertainty coming from the unknown precise value of k 
from NB distribution flts [20] to the pp multiplicity data. In Fig. 3 the forward 
and backward intervals are chosen as: B = {—4,—Ari/2) and F = {Ari/2,4) with 
At] = 0,1,2,3,4,5,6. In Fig. 4 the forward and backward intervals of constant 
widths of 1 are: B = {-Arj/2 - 1, -Ar]/2) and F = {Arj/2,Arj/2 + 1). The pa- 
rameters p and n are calculated using Eqs. (10), (7) and the wounded nucleon 
fragmentation function shown in Fig. 1. All parameters used in these calculations 
are fisted in Tabs. 1 and 2. In both cases the main source of uncertainties is the 
NB parameter k, which is not precisely known for all intervals [20]. 
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Figure 3: The forward-backward multiplicity correlation coefficient 6 as a function 
of the distance A77 between the backward B and forward F intervals. Data points 
(squares) measured in pp at ^/s — 200 GcV arc compared with the results of two 
models: two independent sources of particles (grey band) and the model with a 
single source (dashed band). The widths of the bands reflect the uncertainty in the 
value of k from NB fits to the pp multiplicity data. 
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Figure 4: The same as in Fig. 3 but now the width of each interval is fixed and 
equals 1. 
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As can be observed the model with two independent sources of particles allows 
to understand the main features of the data. It is worth noticing that the strong 
suppression of the correlation coefficient b with increasing Arj is fully determined 
by the suppression of particle production from a single wounded nucleon to the 
backward hemisphere. Clearly, the model in which particles are produced from the 
single source is incorrect. 



Ar] 


F interval 


71 


k 


P 





0.0 < 7] < 4.0 


17.4 


3.70 ±0.30 


0.70 


1 


0.5 < 7] < 4.0 


15.1 


3.85 ±0.35 


0.73 


2 


1.0 < 7] < 4.0 


12.6 


3.95 ±0.40 


0.76 


3 


1.5 <1] < 4.0 


10.2 


4.10 ±0.50 


0.80 


4 


2.0 < 7] < 4.0 


7.71 


4.25 ±0.55 


0.84 


5 


2.5 < 77 < 4.0 


5.39 


4.35 ± 0.60 


0.89 


6 


3.0 < 77 < 4.0 


3.30 


4.50 ±0.70 


0.94 



Table 1: The parameters used in the calculations of the results presented in Fig. 3. 



A7] 


F interval 


71 


k 


P 





0.0 < 7] < 1.0 


4.76 


2.30 ±0.30 


0.54 


1 


0.5 < 7] < 1.5 


4.88 


2.70 ±0.40 


0.59 


2 


1.0 < 77 < 2.0 


4.94 


3.05 ± 0.45 


0.64 


3 


1.5 < 77 < 2.5 


4.78 


3.40 ±0.50 


0.70 


4 


2.0 < 77 < 3.0 


4.41 


3.80 ±0.60 


0.76 


5 


2.5 < 77 < 3.5 


3.90 


4.10 ±0.60 


0.85 


6 


3.0 < 77 < 4.0 


3.30 


4.50 ±0.70 


0.94 



Table 2: The parameters used in the calculations of the results presented in Fig. 4. 

In Fig. 5 the relation (713) \np between the average number of particles (n^) in B 
interval at a given number of particles 7ip in F interval is shown. The measurement 
was performed in two symmetric pseudorapidity intervals B = (—4, 0) and F = 
(0,4). Taking Eqs. (10), (7) into account we obtain p = 0.7 and n = 17.4. The 
main source of uncertainties is the measured value of A; = 3.7 ± 0.3 [20]. 

The model with two independent sources of particles again correctly describes 
the data.^ It is interesting to note that our formalism predicts some deviations from 
linearity, which are too small to be noticeable at a given experimental precision. 

^Except maybe the region of tif ^ 3. 
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Figure 5: The relation between the average number of particles in the backward 
interval {ub) at a given number of particles np va. the forward one. The UA5 pp 
experimental data (squares) at ^/s — 200 GeV are compared with the results of the 
model with two independent sources of particles (grey band) and the model with a 
single source (dashed band). The widths of the bands reflect the uncertainty in the 
value of k from NB fits to the pp multiplicity data. 

5 Conclusions and comments 

Our conclusions can be formulated as follows. 

(i) Assuming that in pp collisions soft particles are produced from two indepen- 
dent sources: left- and right-moving wounded nucleons, we have derived the formulae 
for the forward-backward multiplicity correlation coefficient h and the functional re- 
lation between the average number of particles in the backward interval {ub) at a 
given number of particles n^? in the forward one. This is compared with the case 
where only one source contributes to particle spectrum. 

(ii) We compared our results with the UA5 pp data at ^/s = 200 GeV. We 
conclude that the model with two independent sources of particles allows to under- 
stand the main features of the forward-backward correlation data. As far as the 
correlation coefficient is concerned, we observed very nice qualitative agreement, 
particularly linear suppression of b with increasing distance between the forward 
and backward intervals. This effect is fully determined by the suppression of the 
particle production from a wounded nucleon to the backward hemisphere. 

(iii) We also successfully described the functional relation of the average number 
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of particles in the backward interval (n^) at a given number np of particles in the 
forward one. It is interesting to note that our formalism predicts some deviations 
from linearity, which are too small to be noticeable at a given experimental precision. 
It would be interesting to study this effect in the future experiments. 

(iv) The model in which the particles are produced from a single source is in a 
clear disagreement with the data. 

Following comments are in order. 

(a) The presented analysis was performed only at ^/s = 200 GeV since for higher 
energies the wounded nucleon fragmentation functions are unknown. Studying the 
correlation data at higher energies should allow to extract these functions. 

(b) Assuming (n^) \np to be in the parabohc form with the quadratic term UpC in 
the range np < 40 and np < 50 we obtained c ^ —0.0045 and —0.0035, respectively. 
It is interesting to note that qualitatively similar tendency was observed in the 
quantum-statistical approach [21]. 
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